3.46 \(\int \frac{1}{(a+b \sinh ^2(c+d x))^2} \, dx\)

Optimal. Leaf size=95 \[ \frac{(2 a-b) \tanh ^{-1}\left (\frac{\sqrt{a-b} \tanh (c+d x)}{\sqrt{a}}\right )}{2 a^{3/2} d (a-b)^{3/2}}-\frac{b \sinh (c+d x) \cosh (c+d x)}{2 a d (a-b) \left (a+b \sinh ^2(c+d x)\right )} \]

[Out]

((2*a - b)*ArcTanh[(Sqrt[a - b]*Tanh[c + d*x])/Sqrt[a]])/(2*a^(3/2)*(a - b)^(3/2)*d) - (b*Cosh[c + d*x]*Sinh[c
+ d*x])/(2*a*(a - b)*d*(a + b*Sinh[c + d*x]^2))

________________________________________________________________________________________

Rubi [A]  time = 0.0694698, antiderivative size = 95, normalized size of antiderivative = 1., number of steps used = 4, number of rules used = 4, integrand size = 14, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.286, Rules used = {3184, 12, 3181, 208} \[ \frac{(2 a-b) \tanh ^{-1}\left (\frac{\sqrt{a-b} \tanh (c+d x)}{\sqrt{a}}\right )}{2 a^{3/2} d (a-b)^{3/2}}-\frac{b \sinh (c+d x) \cosh (c+d x)}{2 a d (a-b) \left (a+b \sinh ^2(c+d x)\right )} \]

Antiderivative was successfully verified.

[In]

Int[(a + b*Sinh[c + d*x]^2)^(-2),x]

[Out]

((2*a - b)*ArcTanh[(Sqrt[a - b]*Tanh[c + d*x])/Sqrt[a]])/(2*a^(3/2)*(a - b)^(3/2)*d) - (b*Cosh[c + d*x]*Sinh[c
+ d*x])/(2*a*(a - b)*d*(a + b*Sinh[c + d*x]^2))

Rule 3184

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]^2)^(p_), x_Symbol] :> -Simp[(b*Cos[e + f*x]*Sin[e + f*x]*(a + b*Sin[
e + f*x]^2)^(p + 1))/(2*a*f*(p + 1)*(a + b)), x] + Dist[1/(2*a*(p + 1)*(a + b)), Int[(a + b*Sin[e + f*x]^2)^(p
+ 1)*Simp[2*a*(p + 1) + b*(2*p + 3) - 2*b*(p + 2)*Sin[e + f*x]^2, x], x], x] /; FreeQ[{a, b, e, f}, x] && NeQ
[a + b, 0] && LtQ[p, -1]

Rule 12

Int[(a_)*(u_), x_Symbol] :> Dist[a, Int[u, x], x] /; FreeQ[a, x] &&  !MatchQ[u, (b_)*(v_) /; FreeQ[b, x]]

Rule 3181

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]^2)^(-1), x_Symbol] :> With[{ff = FreeFactors[Tan[e + f*x], x]}, Dist
[ff/f, Subst[Int[1/(a + (a + b)*ff^2*x^2), x], x, Tan[e + f*x]/ff], x]] /; FreeQ[{a, b, e, f}, x]

Rule 208

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-(a/b), 2]*ArcTanh[x/Rt[-(a/b), 2]])/a, x] /; FreeQ[{a,
b}, x] && NegQ[a/b]

Rubi steps

\begin{align*} \int \frac{1}{\left (a+b \sinh ^2(c+d x)\right )^2} \, dx &=-\frac{b \cosh (c+d x) \sinh (c+d x)}{2 a (a-b) d \left (a+b \sinh ^2(c+d x)\right )}-\frac{\int \frac{-2 a+b}{a+b \sinh ^2(c+d x)} \, dx}{2 a (a-b)}\\ &=-\frac{b \cosh (c+d x) \sinh (c+d x)}{2 a (a-b) d \left (a+b \sinh ^2(c+d x)\right )}+\frac{(2 a-b) \int \frac{1}{a+b \sinh ^2(c+d x)} \, dx}{2 a (a-b)}\\ &=-\frac{b \cosh (c+d x) \sinh (c+d x)}{2 a (a-b) d \left (a+b \sinh ^2(c+d x)\right )}+\frac{(2 a-b) \operatorname{Subst}\left (\int \frac{1}{a-(a-b) x^2} \, dx,x,\tanh (c+d x)\right )}{2 a (a-b) d}\\ &=\frac{(2 a-b) \tanh ^{-1}\left (\frac{\sqrt{a-b} \tanh (c+d x)}{\sqrt{a}}\right )}{2 a^{3/2} (a-b)^{3/2} d}-\frac{b \cosh (c+d x) \sinh (c+d x)}{2 a (a-b) d \left (a+b \sinh ^2(c+d x)\right )}\\ \end{align*}

Mathematica [A]  time = 0.302547, size = 96, normalized size = 1.01 \[ \frac{(2 a-b) \tanh ^{-1}\left (\frac{\sqrt{a-b} \tanh (c+d x)}{\sqrt{a}}\right )}{2 a^{3/2} d (a-b)^{3/2}}-\frac{b \sinh (2 (c+d x))}{2 a d (a-b) (2 a+b \cosh (2 (c+d x))-b)} \]

Antiderivative was successfully verified.

[In]

Integrate[(a + b*Sinh[c + d*x]^2)^(-2),x]

[Out]

((2*a - b)*ArcTanh[(Sqrt[a - b]*Tanh[c + d*x])/Sqrt[a]])/(2*a^(3/2)*(a - b)^(3/2)*d) - (b*Sinh[2*(c + d*x)])/(
2*a*(a - b)*d*(2*a - b + b*Cosh[2*(c + d*x)]))

________________________________________________________________________________________

Maple [B]  time = 0.046, size = 749, normalized size = 7.9 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(1/(a+b*sinh(d*x+c)^2)^2,x)

[Out]

-1/d/(tanh(1/2*d*x+1/2*c)^4*a-2*tanh(1/2*d*x+1/2*c)^2*a+4*tanh(1/2*d*x+1/2*c)^2*b+a)*b/a/(a-b)*tanh(1/2*d*x+1/
2*c)^3-1/d/(tanh(1/2*d*x+1/2*c)^4*a-2*tanh(1/2*d*x+1/2*c)^2*a+4*tanh(1/2*d*x+1/2*c)^2*b+a)*b/a/(a-b)*tanh(1/2*
d*x+1/2*c)+1/d/(a-b)/((2*(-b*(a-b))^(1/2)+a-2*b)*a)^(1/2)*arctanh(a*tanh(1/2*d*x+1/2*c)/((2*(-b*(a-b))^(1/2)+a
-2*b)*a)^(1/2))-1/d/(a-b)/(-b*(a-b))^(1/2)/((2*(-b*(a-b))^(1/2)+a-2*b)*a)^(1/2)*arctanh(a*tanh(1/2*d*x+1/2*c)/
((2*(-b*(a-b))^(1/2)+a-2*b)*a)^(1/2))*b-1/d/(a-b)/((2*(-b*(a-b))^(1/2)-a+2*b)*a)^(1/2)*arctan(a*tanh(1/2*d*x+1
/2*c)/((2*(-b*(a-b))^(1/2)-a+2*b)*a)^(1/2))-1/d/(a-b)/(-b*(a-b))^(1/2)/((2*(-b*(a-b))^(1/2)-a+2*b)*a)^(1/2)*ar
ctan(a*tanh(1/2*d*x+1/2*c)/((2*(-b*(a-b))^(1/2)-a+2*b)*a)^(1/2))*b-1/2/d/(a-b)*b/a/((2*(-b*(a-b))^(1/2)+a-2*b)
*a)^(1/2)*arctanh(a*tanh(1/2*d*x+1/2*c)/((2*(-b*(a-b))^(1/2)+a-2*b)*a)^(1/2))+1/2/d/(a-b)/a/(-b*(a-b))^(1/2)/(
(2*(-b*(a-b))^(1/2)+a-2*b)*a)^(1/2)*arctanh(a*tanh(1/2*d*x+1/2*c)/((2*(-b*(a-b))^(1/2)+a-2*b)*a)^(1/2))*b^2+1/
2/d/(a-b)*b/a/((2*(-b*(a-b))^(1/2)-a+2*b)*a)^(1/2)*arctan(a*tanh(1/2*d*x+1/2*c)/((2*(-b*(a-b))^(1/2)-a+2*b)*a)
^(1/2))+1/2/d/(a-b)/a/(-b*(a-b))^(1/2)/((2*(-b*(a-b))^(1/2)-a+2*b)*a)^(1/2)*arctan(a*tanh(1/2*d*x+1/2*c)/((2*(
-b*(a-b))^(1/2)-a+2*b)*a)^(1/2))*b^2

________________________________________________________________________________________

Maxima [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: ValueError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*sinh(d*x+c)^2)^2,x, algorithm="maxima")

[Out]

Exception raised: ValueError

________________________________________________________________________________________

Fricas [B]  time = 2.43671, size = 3767, normalized size = 39.65 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*sinh(d*x+c)^2)^2,x, algorithm="fricas")

[Out]

[1/4*(4*a^2*b - 4*a*b^2 + 4*(2*a^3 - 3*a^2*b + a*b^2)*cosh(d*x + c)^2 + 8*(2*a^3 - 3*a^2*b + a*b^2)*cosh(d*x +
c)*sinh(d*x + c) + 4*(2*a^3 - 3*a^2*b + a*b^2)*sinh(d*x + c)^2 + ((2*a*b - b^2)*cosh(d*x + c)^4 + 4*(2*a*b -
b^2)*cosh(d*x + c)*sinh(d*x + c)^3 + (2*a*b - b^2)*sinh(d*x + c)^4 + 2*(4*a^2 - 4*a*b + b^2)*cosh(d*x + c)^2 +
2*(3*(2*a*b - b^2)*cosh(d*x + c)^2 + 4*a^2 - 4*a*b + b^2)*sinh(d*x + c)^2 + 2*a*b - b^2 + 4*((2*a*b - b^2)*co
sh(d*x + c)^3 + (4*a^2 - 4*a*b + b^2)*cosh(d*x + c))*sinh(d*x + c))*sqrt(a^2 - a*b)*log((b^2*cosh(d*x + c)^4 +
4*b^2*cosh(d*x + c)*sinh(d*x + c)^3 + b^2*sinh(d*x + c)^4 + 2*(2*a*b - b^2)*cosh(d*x + c)^2 + 2*(3*b^2*cosh(d
*x + c)^2 + 2*a*b - b^2)*sinh(d*x + c)^2 + 8*a^2 - 8*a*b + b^2 + 4*(b^2*cosh(d*x + c)^3 + (2*a*b - b^2)*cosh(d
*x + c))*sinh(d*x + c) - 4*(b*cosh(d*x + c)^2 + 2*b*cosh(d*x + c)*sinh(d*x + c) + b*sinh(d*x + c)^2 + 2*a - b)
*sqrt(a^2 - a*b))/(b*cosh(d*x + c)^4 + 4*b*cosh(d*x + c)*sinh(d*x + c)^3 + b*sinh(d*x + c)^4 + 2*(2*a - b)*cos
h(d*x + c)^2 + 2*(3*b*cosh(d*x + c)^2 + 2*a - b)*sinh(d*x + c)^2 + 4*(b*cosh(d*x + c)^3 + (2*a - b)*cosh(d*x +
c))*sinh(d*x + c) + b)))/((a^4*b - 2*a^3*b^2 + a^2*b^3)*d*cosh(d*x + c)^4 + 4*(a^4*b - 2*a^3*b^2 + a^2*b^3)*d
*cosh(d*x + c)*sinh(d*x + c)^3 + (a^4*b - 2*a^3*b^2 + a^2*b^3)*d*sinh(d*x + c)^4 + 2*(2*a^5 - 5*a^4*b + 4*a^3*
b^2 - a^2*b^3)*d*cosh(d*x + c)^2 + 2*(3*(a^4*b - 2*a^3*b^2 + a^2*b^3)*d*cosh(d*x + c)^2 + (2*a^5 - 5*a^4*b + 4
*a^3*b^2 - a^2*b^3)*d)*sinh(d*x + c)^2 + (a^4*b - 2*a^3*b^2 + a^2*b^3)*d + 4*((a^4*b - 2*a^3*b^2 + a^2*b^3)*d*
cosh(d*x + c)^3 + (2*a^5 - 5*a^4*b + 4*a^3*b^2 - a^2*b^3)*d*cosh(d*x + c))*sinh(d*x + c)), 1/2*(2*a^2*b - 2*a*
b^2 + 2*(2*a^3 - 3*a^2*b + a*b^2)*cosh(d*x + c)^2 + 4*(2*a^3 - 3*a^2*b + a*b^2)*cosh(d*x + c)*sinh(d*x + c) +
2*(2*a^3 - 3*a^2*b + a*b^2)*sinh(d*x + c)^2 - ((2*a*b - b^2)*cosh(d*x + c)^4 + 4*(2*a*b - b^2)*cosh(d*x + c)*s
inh(d*x + c)^3 + (2*a*b - b^2)*sinh(d*x + c)^4 + 2*(4*a^2 - 4*a*b + b^2)*cosh(d*x + c)^2 + 2*(3*(2*a*b - b^2)*
cosh(d*x + c)^2 + 4*a^2 - 4*a*b + b^2)*sinh(d*x + c)^2 + 2*a*b - b^2 + 4*((2*a*b - b^2)*cosh(d*x + c)^3 + (4*a
^2 - 4*a*b + b^2)*cosh(d*x + c))*sinh(d*x + c))*sqrt(-a^2 + a*b)*arctan(-1/2*(b*cosh(d*x + c)^2 + 2*b*cosh(d*x
+ c)*sinh(d*x + c) + b*sinh(d*x + c)^2 + 2*a - b)*sqrt(-a^2 + a*b)/(a^2 - a*b)))/((a^4*b - 2*a^3*b^2 + a^2*b^
3)*d*cosh(d*x + c)^4 + 4*(a^4*b - 2*a^3*b^2 + a^2*b^3)*d*cosh(d*x + c)*sinh(d*x + c)^3 + (a^4*b - 2*a^3*b^2 +
a^2*b^3)*d*sinh(d*x + c)^4 + 2*(2*a^5 - 5*a^4*b + 4*a^3*b^2 - a^2*b^3)*d*cosh(d*x + c)^2 + 2*(3*(a^4*b - 2*a^3
*b^2 + a^2*b^3)*d*cosh(d*x + c)^2 + (2*a^5 - 5*a^4*b + 4*a^3*b^2 - a^2*b^3)*d)*sinh(d*x + c)^2 + (a^4*b - 2*a^
3*b^2 + a^2*b^3)*d + 4*((a^4*b - 2*a^3*b^2 + a^2*b^3)*d*cosh(d*x + c)^3 + (2*a^5 - 5*a^4*b + 4*a^3*b^2 - a^2*b
^3)*d*cosh(d*x + c))*sinh(d*x + c))]

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*sinh(d*x+c)**2)**2,x)

[Out]

Timed out

________________________________________________________________________________________

Giac [A]  time = 1.21252, size = 196, normalized size = 2.06 \begin{align*} \frac{{\left (2 \, a - b\right )} \arctan \left (\frac{b e^{\left (2 \, d x + 2 \, c\right )} + 2 \, a - b}{2 \, \sqrt{-a^{2} + a b}}\right )}{2 \,{\left (a^{2} d - a b d\right )} \sqrt{-a^{2} + a b}} + \frac{2 \, a e^{\left (2 \, d x + 2 \, c\right )} - b e^{\left (2 \, d x + 2 \, c\right )} + b}{{\left (a^{2} d - a b d\right )}{\left (b e^{\left (4 \, d x + 4 \, c\right )} + 4 \, a e^{\left (2 \, d x + 2 \, c\right )} - 2 \, b e^{\left (2 \, d x + 2 \, c\right )} + b\right )}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*sinh(d*x+c)^2)^2,x, algorithm="giac")

[Out]

1/2*(2*a - b)*arctan(1/2*(b*e^(2*d*x + 2*c) + 2*a - b)/sqrt(-a^2 + a*b))/((a^2*d - a*b*d)*sqrt(-a^2 + a*b)) +
(2*a*e^(2*d*x + 2*c) - b*e^(2*d*x + 2*c) + b)/((a^2*d - a*b*d)*(b*e^(4*d*x + 4*c) + 4*a*e^(2*d*x + 2*c) - 2*b*
e^(2*d*x + 2*c) + b))